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Abstract 

In many advanced low NO x gas turbine combustion techniques, such as rich-burn/quick- 
mix/lean-burn (RQL), jet mixing in a reacting, hot, fuel-rich crossflow plays an important role in 
minimizing all pollutant emissions and maximizing combustion efficiency. Assessing the degree 
of mixing and predicting jet penetration is critical to the optimization of the jet injection design 
strategy. Different passive scalar quantities, including carbon, oxygen, and helium are compared 
to quantify mixing in an atmospheric RQL combustion rig under reacting conditions. The results 
show that the 0 2 -based jet mixture fraction underpredicts the C-based mixture fraction due to jet 
dilution and combustion, whereas the He tracer overpredicts it possibly due to differences in 
density and diffusivity. The He-method also exhibits significant scatter in the mixture fraction 
data that can most likely be attributed to differences in gas density and turbulent diffusivity. The 
jet mixture fraction data were used to evaluate planar spatial unmixedness, which showed good 
agreement for all three scalars. This investigation suggests that, with further technique 
refinement, 0 2 or a He tracer could be used instead of C to determine the extent of reaction and 
mixing in an RQL combustor. 


Nomenclature 

d orifice axial length 

/ jet mixture fraction 

f Y i jet mixture fraction based on mass fraction of species i 

fxi jet mixture fraction based on molar fraction (or concentration) of species i 

f var variance of all point / values in a plane w.r.t .f avg 

f avg area- weighted average jet mixture fraction specific to each plane 

J jet to crossflow momentum-flux ratio 

Mi molar mass of species i 

§ equivalence ratio = (fuel/air)i 0ca i/(fuel/air) st0 i ch i 0metric 


Corresponding author. 
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R radius of the quick-mix module 

U s spatial unmixedness 

Xi molar fraction of species i 

Yt mass fraction of species i 

x axial distance from the leading edge of the orifices 


Introduction 

Many advanced low NO x combustion techniques, such as lean premixed prevaporized 
injection, lean direct injection, and rich-burn/quick-mix/lean-bum (RQL) rely on the rapid and 
thorough mixing of air and fuel in order to minimize all pollutant emissions and maximize 
combustion efficiency. Various studies have found that fuel-air ratio non-uniformities 
significantly affect NO x emissions [1,2,3]. 

In gas turbine combustion, jet mixing in a reacting, hot, fuel-rich crossflow plays an important 
role due to air jet injection in the primary, secondary, and dilution zones of the combustor. 
Assessing the degree of mixing and predicting jet penetration is especially critical in the RQL 
combustion concept. One of the advantages of RQL over other combustion techniques is quick 
and complete mixing between the rich and lean zones of the combustor in order to eliminate hot, 
near-stoichiometric reactant pockets that may lead to NO x formation. In addition to combustion, 
the assessment of jet mixing into a crossflow can be applied to a wide range of fields such as gas 
turbine cooling and staging, fuel-air premixing, vertical short takeoff and landing aircraft, and 
pollutant discharge from stacks or pipes. 

This study compares the use of different scalar quantities, including carbon, oxygen, and an 
inert tracer gas to quantify mixing in an atmospheric RQL combustion rig under reacting 
conditions. 


Background 

Most experimental jet-in-crossflow studies have focused on non-reacting systems, with only a 
limited number of tests having been reported under reacting conditions. Although isothermal 
testing is useful and convenient, actual combustor mixing and performance needs to be measured 
in a combusting flow. An extensive listing of these isothermal and reacting studies can be found 
in [4, 5, 6, 7]. 

The diagnostic technique chosen to determine mixing in reacting systems is important to the 
outcome of this study. Qualitative characterization of the mixing process can be inferred through 
the measurement of temperature profiles and species concentrations. However, to determine the 
true extent of mixing in a reactant flow field, one needs to measure, experimentally, the jet 
mixture fraction f 

Two techniques, used in numerous studies and summarized by Jones et al. [8], can be used to 
measure or deduce / in a transverse flow. The first method is to use non-intrusive optical 
diagnostics, such as laser-induced fluorescence (LIF), Rayleigh scattering [9,10], or Raman 
scattering, to quantify spatial and temporal unmixedness. Planar imaging eliminates the need for 
extractively measuring multiple species and for sampling at multiple points in the flame. 


NAS A/CR— 2004-2 12886 


2 



However, optical methods are limited by the availability of costly laser and detector systems, 
the practicality of optical access into the flow, the need to seed or modify the fuel to obtain 
the desired optical signal, as well as in-flame interferences such as molecular quenching. 

The second approach, which is the focus of this study, is based on the measurement of a 
conserved scalar (i.e., unaffected by the chemical reaction), such as carbon, oxygen, nitrogen, 
or hydrogen mass fraction, equivalence ratio, or an inert gas. This treatment assumes that the 
slowest chemical kinetic reaction rate is much faster than the turbulent mixing time scales [11]. 
Under this assumption, local instantaneous composition measurements correspond to chemical 
equilibrium and can be related to a strictly conserved scalar variable. Furthermore, this technique 
assumes a well-mixed recirculation zone in which mixing times are much less than residence 
times [11]. 

Previous studies [8,12,13] have used carbon mass fraction and equivalence ratio to calculate f 
Aspirated emissions samples were analyzed for CO, C0 2 , 0 2 , total hydrocarbons (THCs), and, 
in the case of Jones et al. [8], H 2 . These analyses directly capture all carbon-carrying species 
needed to determine a C-based f However, in the case of a jet in a rich crossflow, other 
quantities, such as the oxygen atom (O), the oxygen molecule (0 2 ), or an inert tracer gas such 
as helium (He) or neon (Ne), can be used to simplify the calculation of f With respect to 0 2 , 
because it is a key participant in the combustion reactions, it cannot truly be considered a 
conserved scalar. However, since it is only present in the jet flow and since there is considerable 
excess air present in the downstream flow, its concentration profile can not only be used to 
indicate jet presence and dispersion, as was shown in [12,13], but potentially mixture fraction 
as well. 

Helium has served as an inert tracer in various studies to determine a wide range of 
parameters, such as groundwater transport [14], fluid flow in a porous rock [15], residence time 
in a spray-drying tower [16], automobile exhaust flow rate [17], impervious wall effectiveness 
of film-cooling slots [18], scalar flow field in a combustor rig under isothermal (non-reacting) 
conditions [1 1,19,20], and mass transport rates in a non-reacting jet-in-crossflow [21,22]. Helium 
is an inexpensive, readily available gas that is detectable using gas chromatography, mass 
spectrometry, or a catharometer (thermal conductivity detector). 

The purpose of this study is to use alternative passive scalars, namely O, 0 2 , and He to 
generate jet mixture fractions at specific planes in an RQL combustion rig. The mixture fraction 
data are then used to determine the degree of spatial mixing of chemical species at each of the 
measurement planes in the combustor. The results are compared to carbon-based / results to 
demonstrate the viability of a simpler method that requires the analysis of only a single 
compound to quantify air-fuel mixing, and to determine suitability of this diagnostic in this 
experiment and for future reacting applications. 

This study expands upon two previous studies by Leong et al [12,13]. In ref. 12 the setup 
for reacting tests in an RQL crossflow configuration was described and characterized, and in 
ref. 13 the optimal number of jet injection orifices was determined to obtain rapid mixing of air 
jets in a rich crossflow and a uniformly lean, low temperature mixture at the exit plane of the 
combustor. Most of the species concentrations required for the carbon and oxygen-based mass 
mixture fractions calculations were collected during the second study [13]. In the work reported 
herein, a protocol for the sampling and analysis of the helium tracer gas is established. Tracer gas 
concentrations are then measured under reacting conditions for a series of RQL modules with a 
different number of jet orifices. Data are collected at specified planes and spatial coordinates to 
allow direct comparison with results obtained in ref. 13. 
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Experiment 


This section describes the experimental setup, the data measurement protocol, and the 
procedure for calculating carbon-based and tracer gas-based mixture fractions. 


Reacting Test Facility 

The reacting jet-in-crossflow experimental set-up, shown in figure 1, has been described in 
detail in previous papers [12,13]. The upward-fired atmospheric test facility supplies a uniform, 
fuel-rich composition of gases to the quick-mix section. The quick-mix section utilizes 
interchangeable quartz tubes containing different jet orifice configurations. 

A schematic of the quartz module, shown with the location of the planes of interest, is 
depicted in figure 2. The inner and outer diameters of the tube are, respectively, 80 mm by 
85 mm, and its length is 280 mm. Four different modules were tested, with each one containing 
a different number of circular orifices (10, 12, 14, or 18) arranged equidistantly around the 
circumference of the tube. (Note that measurement results for an 8 hole module were reported in 
ref. 13, but are not included here because He data were not collected for this module). The orifice 
centerlines are located 115 mm from the entrance of the module. The four measurement planes 
are displaced from the orifice leading edge as follows: one duct radial length upstream ( x/R = -1), 
one orifice diameter downstream ( x/R = d/R ), and one-half and one duct radial lengths 
downstream (x/R= 0.5, 1). This particular region was chosen because, based on data presented 
by Leong et al. [13], this is where the greatest changes in spatial unmixedness occurred. For 
each module tested, point samples were taken from a two-orifice sector at 16 radially equidistant 
locations (see figure 2). 

Emissions data, required for calculation of carbon and oxygen-based mixture fractions, were 
previously presented in ref. 13 for planes x/R = -1, d/R , and 1. Additional emissions samples, 
using the same setup and procedure as outlined in refs. 12 and 13, were collected at x/R = 0.5 
because intense mixing activity occurs at this plane according to tracer gas data. Emissions were 
analyzed using the following techniques: CO and C0 2 by non-dispersive infrared absorption, 0 2 
by paramagnetism, and THCs by flame ionization. 

The experiment utilizes gaseous propane as the combustion fuel. Propane is first mixed with 
air to yield a fuel-air equivalence ratio $ of 1 .67 in the rich combustion section. The rich product 
generation is described in more detail in ref. 12. The rich-burning mixture, with average 
temperatures at the x/R= -1 plane of 1500 K, enters the quartz mixing section and undergoes 
additional reaction with jets of air to result in an overall $ of 0.45. The jet air is fed by a plenum 
that surrounds the rich combustion chamber and the quartz tube. Heat transferred from the 
combustor to the plenum air heats the jets of air to 480 K before they enter the jet-mixing section. 

The rich equivalence ratio is obtained by setting the mass flow rates of propane and crossflow 
air at 2.96 g/sec and 27.5 g/sec, respectively. The lean equivalence ratio downstream of the 
jets requires a total jet mass flow rate of 75.2 g/sec. The reference velocity of the total flow is 
18 m/sec. Based on the temperatures measured in the reacting system, the jet-to-crossflow 
density ratio is 3.3, and the jet-to-crossflow momentum-flux ratio ./is 57. This set of conditions 
is the same as utilized in previous experiments [12,13], and was selected to fall within the range 
of gas turbine combustor operating conditions [4,5,6,23]. Jis kept constant by keeping the total 
effective orifice area (903 mm 2 ) constant for each of the modules tested. This results in different 
orifice diameters for each orifice number configuration, namely 12.5, 1 1.5, 10.6, and 9.4 mm 
respectively for the 10, 12, 14, and 18 round hole modules. 
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The Mixture Fraction Determination 


The experiment consists of jet air streams injected into a cylindrically-confined crossflow of a 
fuel-rich mixture of partially-reacted propane and air. Based on the method by Jones et al. [8], 
the mixture fraction in terms of the jet fluid is defined as 


y crossflow y sample 

f =— — 

J i Y cross fl° w y J et 


(i) 


where the jet mixture fraction /is related in terms of the mass fractions /• of a conserved scalar i 
in the crossflow, jets, and extracted gas sample. The jet mixture fraction tracks the amount of 
the jet fluid relative to the total mixture at a specific location in the combustor. A gas sample 
composed entirely of crossflow fluid yields a value of f= 0, while a sample comprised wholly 
of jet fluid produces a value of f= 1. 

As shown in table 1, four different formulations for mass-based jet mixture fraction were 
examined. In each case, the table lists the measured quantities, the minimum number of unknown 
variables needed to calculate / and the assumptions and equations required to solve for these 
unknowns. The mass-fraction based calculations using C, O, 0 2 , and He as the conserved scalars 
follow the procedure outlined in [8]. Key assumptions in this method are: 

1. The combustion gas mixture is composed of major species only (i.e., CO, C0 2 , 0 2 , 

H 2 0, N 2 , H 2 , and THCs) and the respective molar fractions sum to unity. 

2. The 0 2 /N 2 and thus O/N ratios are the same in the sample stream as in the combustion air 

3. The sampled C/H molar ratio is the same as in the fuel stream. 

4. To ensure a unique solution to the system of linear equations (i.e., equal number of 
equations and unknowns), two hydrocarbons are included: C 3 H 8 as well as C 2 H 4 , which is 
a byproduct from the pyrolysis of C 3 H 8 . 

5. In the absence of a direct measurement for hydrogen, the H 2 molar fraction is assumed to 
be proportional to the CO molar fraction [24,25]. 

Assumptions 2 and 3 also imply that all major species have equal diffusivities, a reasonable 
assumption in turbulent flows according to [1 1]. Assumption 3 is used to infer a molar fraction X 
for H 2 0 and thus convert the measured dry basis emissions X i dry to a wet basis X Uwet to represent 
the gases as found in the actual combustion reaction. To perform this basis conversion, Jones et 
al. [8] uses a matrix formulation of the form 


v4X wet X d 


ry 


with a solution given by 


x wet =A- 1 x l 


dry 


( 2 ) 

( 3 ) 


where A is an invertible, square matrix, representing the coefficient matrix in the linear system of 
equations setup using the assumptions listed in table 1. 

The C- and O-based mixture fractions are considered to give the most comprehensive / values, 
since the calculations use all available data and only the five assumptions just listed. In fact, 
the values for f YC and f Y o will be almost equal given that both fractions are calculated using the 
same linear equations set under the same assumptions. Thus, all f YC results in this study can be 
considered interchangeable with f m . It should be noted that oxygen atom-based / calculations 
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may be subject to higher percentage uncertainties under certain assumptions due to the 
subtraction of two nearly equal numbers (i.e., in the numerator of the mixture fraction formula) 
[26]. 

Using 0 2 and He as passive scalars can simplify the determination of f since only data from a 
single species is needed for the calculation if molar fractions are used. In the case of 0 2 , under 
equilibrated lean combustion conditions, the concentration will vary as a direct function of 
equivalence ratio and, by definition of f of mixture fraction as well. Any 0 2 present in rich 
regions can be considered as a diluent, which also allows it to serve as quasi-conserved scalar for 
determining mixture fraction. As shown in table 1, the ratio of measured to maximum (20.9%) 

0 2 concentrations can be used to represent 

In the tracer gas case, the helium was injected into the crossflow air rather than into the jets 
to maintain a fair comparison between the mixing fields obtained by the previous carbon atom 
tracking and the inert tracer method. (In the carbon-based method, the only source of carbon is 
from the crossflow, not the air jets). Seeding the jet flow, which is 2.5 times the mass flow of 
the crossflow, would have required helium flow rates that would have quickly depleted available 
helium supplies during the course of a test. In addition, seeding the crossflow allows for the 
assessment of its uniformity. 

Initially, neon and argon were considered as tracers, since their molecular weights and thus 
densities and diffusivities more closely match those of propane and air as shown in table 2. 

Argon, however has a very high background concentration (9340 ppm) compared to helium 
(5 ppm) and neon (18 ppm) [27] and thus was not pursued. Preliminary studies were conducted 
on the RQL module using the procedure in order to compare helium and neon as inert tracers. 
Concentration measurements for both gases injected into the crossflow yielded similar results, 
suggesting that either could be used as a tracer. Helium has an extensive history as an inert tracer 
and is much cheaper and more readily available than neon. These practical considerations lead to 
the choice of helium for use in this experiment. 

The injection of the helium tracer into the crossflow, the assumption that the jet air stream 
contains negligible levels of helium, and the approximation that the helium concentration at the 
crossflow injection plane is relatively uniform lead to the helium-based jet mixture fraction 
relationship shown in table 1. (Variation across the plane was found to be less than 3%.) 

The value f X He gives the relative change in helium concentration with respect to the maximum 
concentration in the plane x/R = -1 . Helium gas of ultra-pure carrier grade (99.999% purity) is 
supplied at a flow rate of 4.5 1/min and injected into the propane stream. This is the minimum 
flow rate tested under the current operating conditions that produces a distinct signal, which 
corresponds to a maximum volume concentration of 0.3% of the rich crossflow mixture. The 
propane -helium mixture is injected into the crossflow air, and flows through a 4.3-m mixing 
length filled with baffles to prepare the gas mixture for combustion. The resulting mass fractions 
of helium, fuel, and air upstream and downstream of the jet mixing section are noted in table 3. 

First, the mass-fraction-based f Y0 2 and f Y He were determined in order to provide a point of 
reference to compare to f YC and f YQ . Then, the molar fraction (or concentration)-based/ra 2 and 
fxHe were calculated. The accuracy of using molar fraction instead of mass fraction relies on the 
fact that most of the sample mass and volume at the measurement planes within the jet injection 
section are composed of jet air. As a result, one can make the following approximations: 

1 . The total mass of the sample is constant (it actually varies by ±4% across the sampling 
planes) 

2. X dry is linearly proportional to X wet (the average difference between X dry and X wet is 
approximately 6% for 0 2 and He). 
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3. 


Because of approximations 1 and 2 and given the linear relationship between ) and X,- 


(i.e., Y, = 


M 


M 


-X ( , where M t is the molar mass of species i and M samp i e is the 


sample 


molecular weight of the sample), fn, W et =fxi,dry • The maximum deviation off x with 
respect to f Y was calculated to be ±0.08 forf xc , ±0.03 for f X 0 2, and ±0.05 for f XHe . 


The procedure for collecting the CO, C0 2 , 0 2 , and THC emissions used to determine f YC ,fYo, 
f Y o2, and f X0 2 was outlined in refs. 12 and 13. The experimental component in this paper involves 
the measurement of the helium tracer gas used for calculating f YHe andf XHe . 


Inert Gas Tracer Sampling Protocol 

The injection and sampling train of the helium tracer system is shown in figure 3. Gas 
samples are extracted from the flow field and drawn through a water-cooled probe with a pump. 
After condensing water from the sample through an impinger submersed in an ice bath, the gas 
is sent to a gas chromatograph (HP 5890 Series II) to measure the helium concentration. When 
a sample analysis is initiated, the volume of gas contained in the 250-pl sample loop of the gas 
chromatograph is injected into the column, while the balance of the extracted gas is diverted 
through a flow bypass. 

The gas chromatograph setup is optimized to obtain the fastest elution time of the helium 
tracer gas with acceptable chromatographic separation. In order to separate the helium atoms 
from the heavier molecules of fuel, air, and combustion products, the sample flows through 
two columns, connected in series. Both columns, manufactured by J&W Scientific, Inc., are 
megabore, capillary, gas solid (GS) columns, with an inner diameter of 0.54 mm, and a length of 
30 m. Ultra-pure grade hydrogen (99.999% H 2 concentration) is used as the carrier gas in order 
to minimize the sample elution time, as well as prevent chromatographic interference from H 2 
generated from the combustion of fuel-rich mixtures. 

Chromatographic separation takes place in two stages. The gas sample is injected into the 
GS-Q column before flowing through the GS-MolSieve column. The GS-Q column separates 
hydrocarbon molecules such as methane, ethane, and propane, while the GS- MolSieve column 
separates compounds of low molecular weight, including helium, neon, oxygen, nitrogen, and 
carbon monoxide. The effects of the hydrocarbon molecules on the sensitivity of the MolSieve 
column necessitated the use of timed valve switching. The valve controlling the flow circuitry is 
switched 30 sec after the helium enters the GS-MolSieve column and interrupts the flow between 
the GS-Q and GS-MolSieve columns. The carrier gas flow through the GS-Q column is reversed 
(backflushed) to flush the sample containing the hydrocarbon molecules out of the system, while 
the helium continues to separate from the remainder of the sample in the GS-MolSieve column as 
it moves toward the detector. 

To detect helium, the gas chromatograph uses a thermal conductivity detector (TCD). The 
TCD detects the difference in the thermal conductivities of the eluted sample and the carrier gas, 
generating a differential voltage signal. The output from the TCD is connected to an integrator 
(Spectra Physics DataJet), which in turn is connected to a computer. The WINner on Windows 
software package by Thermo Separation Products was used to integrate the resulting 
chromatograms. From helium calibration runs performed before and after each test, a constant of 
proportionality is obtained to quantify the integrated areas under the peaks in terms of volumetric 
concentration. 

Completion of the analysis requires approximately 1.5 min, followed by 30 sec to flush out the 
columns to prepare them for the next sample. The gases elute in the following order: helium at 
0.59 min, air at 0.66 min, and CO at 0.74 min. 
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Results and Discussion 


Mixture Fraction 

The simplified molar mixture fraction formulations proposed in this paper, namely f X o 2 and 
fxHo are plotted in figure 4 and 5 versus f YC and equivalence ratio (|). Figure 4 also graphs^ 
versus f YC , showing that the molar fraction-based mixture fraction slightly overpredicts f YC by 
approximately 5%. The equivalence ratio was calculated from f YC according to Jones et al. [11] 
using the following formula: 


where 


9 = 


4 




i+- 


Y n 


5 Mr 


M. 


} 2 J ~ 


4 = 


Y sample 


Y, 


fuel 


= (1 -/) 


M, 




M 


crossflow 


(4) 


(5) 


Figures 4 and 5 include data collected at all the measurement planes for the four test modules. 
The data are grouped in terms of sampling radius to better explain certain results. Grouping by 
module type or sampling plane did not yield any identifiable trends. Figures 4 and 5 show that 
fxo 2 and f X He correlate positively with / 7C , but they exhibit significant deviations from the carbon- 
based mixture fraction. Experimental uncertainties for fxHe,fxo 2 , and f YC are, respectively, ±0.08, 
±0.01, and ±0.02. The apparent gap in data points in the region of 0.2 </< 0.4, which 
corresponds to near-stoichiometric combustion, illustrates the staging in RQL combustion that 
causes the reaction to shift from § = 1.67 to § = 0.45 by quick mixing to avoid hot stoichiometric 
pockets. 

In the case of fxo 2 , the values underpredict f YC . Two trends are evident in figure 4, one 
corresponding to combustion equilibrium (predicted by the NASA Equilibrium Code) and the 
other to jet dilution and quenching of the reaction. The dilution effect is linear in both the slightly 
rich and slightly lean regimes. Extra scatter in the data around the knee at near-stoichiometric 
conditions (i.e., for 0.3 </< 0.5) indicates non-equilibrium conditions, resulting in higher than 
equilibrium 0 2 and CO values. Figure 4 suggests that the deviation occurs away from the walls 
(i.e., r = R/ 3, 27?/3), from which one can surmise that the non-equilibrium conditions are probably 
due to the intense, turbulent, quick mix transition from rich to lean conditions following jet 
injection. Near the edges of the combustor (i.e., r = R ), where jet interaction is less, equilibrium 
conditions prevail. Thus, in addition to the usefulness of 0 2 concentration profiles to track jet 
trajectories and determine the general zones where reacting and mixing processes are occurring 
(as shown in ref. 12), one can use 0 2 to distinguish between the two processes. 

As seen in figure 5, thzfxHe data overpredict and show considerable scatter with respect to f YC . 
The overprediction may be attributed to (a) the higher molecular and turbulent diffusivity of He 
with respect to the other major products of combustion, and (b) the use of molar fraction rather 
than mass fraction which can introduce a positive bias of up to 0.05. for high mixture fraction 
values. The data collected at r = R, where, as previously noted, equilibrium conditions hold, 
appear to correlate well with f YC . However, once the jet interaction is thrust into the picture, the 
tracer may be unable to follow the high momentum jet trajectories. In ref. 12, it was shown that 
overpenetrating jets displace the rich reacting fluid towards the walls, while the jet mass migrates 
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and accumulates in the central core of the combustor. This is corroborated by the fact that the 
flows at r = 2R/3 and R / 3 seem to either completely trap (f= 0) or totally exclude (f= 1) the He 
tracer molecules. 

The scatter, on the other hand, may be attributed to a variety of factors, in addition to 
experimental uncertainty. The high turbulence at the jet-crossflow interface due to reacting and 
mixing processes can cause large spatial-temporal fluctuations and asymmetries in the mixture 
fraction. Density differences between He and the other major combustion products can result in 
large differences in momentum flows and kinematic eddy viscosities (turbulent momentum 
diffusivities), which in turn affects the He transport in the turbulent eddies in the mixing regions. 
Previous isothermal experiments used He concentrations of at least 1%, while, in this study, the 
concentrations were less than 0.3%. Although these values were chosen to minimize a potential 
diluent effect, perhaps, in hindsight, the concentration levels were too low for use in reacting 
conditions. Finally, there is the possibility of small offsets in the relative position of the sampling 
probe between runs. As a result of the extreme scattering in the He data, this technique would 
need further refinement before it would be seen as amenable for reliably determining mixture 
fraction. 

To better quantify linear correlations between the mixture fraction data, two parameters 
were selected: (1) standard errors (a), which give the magnitude of typical deviation from the 
estimated linear regression line, and (2) coefficients of determination (R 2 ), which represent the 
proportion of a data set that can be explained by the linear regression model. These values are 
presented in table 4. Thus, a and R 2 quantify the data scatter, while the regression curve fit 
would quantify the actual agreement between the various mixture fraction formulations. The 
statistical analysis confirms the better fit for the 0 2 data, whose a is on the order of half of that of 
the He data. The analysis also shows the He data to be consistently overpredicting/ 7C by -10% 
and 0 2 underpredicting f YC also by -10%. 

To better visualize and compare the various scalars, stacked sector plots, presented in figure 6 
to 8, were generated to show the evolution of the scalar flow field through the modules in terms 
of f YC (reproduced from ref. 13), f X o 2 , and f X He- Note that the flow is upward such that the farthest 
upstream plane is at the bottom, and the farthest downstream plane is at the top. The f X o 2 and f XHe 
results suggest a qualitative match with the f Y c data. All plots show uniformity at x/R = -1 and 
similar trends in jet penetration at x/R = d/R. 

The ideal, fully-mixed case occurs when the jets mix uniformly with the crossflow. One can 
calculate the ideal jet mixture fraction by substituting the jet and crossflow mass fractions, based 
on input mass flow rates, into eq. (1) to yield an / of 0.714. The f Y c stack plots in figure 6 indicate 
that such a region occurs at the plane x/R = 1. Th ef X He data shown in figure 8 indicate poorer 
mixing with most of the downstream flow composed of jet fluid. The f X o 2 data plotted in figure 7, 
on the other hand, suggest only a narrow band of good mixing, but with crossflow fluid still 
present in the central and outer rings of the combustor. Both helium and oxygen-based / data 
show steeper gradients than their equivalent carbon-based values. 


Spatial Unmixedness 


To determine the effectiveness of the tracer gas method in assessing overall mixing, one can 
calculate spatial unmixedness U s based on f U s is the normalized variance quantifying planar 
mixing and is defined by as: 


U s = 



( 6 ) 
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where f var refers to the variance of all/s in a plane that deviate from f avg , the area- weighted 
average jet mixture fraction specific to each plane [28]. U s values lie between 0 (perfect fuel- 
air mixing) and 1 (totally unmixed system). Figure 9 reveals that U s ,xo2 and U s ,xHe correlate 
reasonably well with each other and with U s ,yc- The same data are plotted versus axial distance 
x/R in figure 10. At the orifice trailing edge (i.e., where x/R = d/R ) where all the jet mass is 
injected, U s reaches a peak and then decreases downstream as the jets mix with the crossflow. 

It appears that the procedure of calculating a planar variance and then normalizing the result 
eliminates some of the scatter and discrepancies seen in the point measurement mixture fraction 
data (see table 4 for the corresponding a and R 2 regression values). 


Comparison of Conserved Scalar Measurement Methodologies 

The quantitative results discussed above warrant a brief discussion of the relative advantages 
and disadvantages of each technique. The carbon-based method requires several emissions 
analyzers for simultaneous species measurement, but also permits localized characterization of 
pollutant emissions, which can then be evaluated with respect to local mixing efficiency. This 
method is also the most comprehensive technique since all major species are included in the 
mixture fraction calculation. 

The oxygen-based method only uses a single (0 2 ) analyzer, but would require a back 
calculation to correct for dilution or lean combustion in order to extract the actual local mixture 
fraction. However, this method appears to be the simplest technique for the rapid characterization 
of local jet mixture fraction and spatial unmixedness, despite the fact that the effective use of 0 2 
as a passive scalar is predicated on the assumption that 0 2 is present in only one of the two flows. 
Both the C-based and 0 2 -based methods allow one to assess the degree of reaction versus mixing 
occurring between the two streams. 

The helium-based method, on the other hand, only provides mixing information and requires 
the added complication of metering and injecting a tracer gas. However, He-based measurements 
could serve as simpler diagnostic substitute to quantify mixing in the absence of a full emissions 
measurement console. Furthermore, this method could be used to independently verify the 
carbon-based mixture fraction calculations. The gas chromatographic analysis used for He 
detection could also be configured to measure emissions and helium simultaneously. Hence, 
using a single sampling system, one could then apply any of the above methods to characterize 
the mixing field. 

If helium or other inert tracer gases are to be used in future experiments to determine mixture 
fractions, the effect of tracer gas density and diffusivity on its dispersion in the mixing field need 
to be investigated. It would be useful to determine the linearity of the tracer injection method, 
since a normalized response curve should be independent of quantity of tracer injected [16]. In 
addition, the sensitivity and repeatability of the tracer gas analysis system should be improved, 
for example, by ensuring a constant pressure sample injection and using multipoint calibration 
for the gas chromatographic analysis. Use of a commercial helium detector, instead of gas 
chromatography, would also help to overcome sensitivity issues, as well as speeding up the 
gas analysis. 


Conclusions 

This study investigated the use of passive scalars, namely the carbon atom, the oxygen 
molecule, and helium (as an inert tracer gas) to quantify jet mixing in a reacting crossflow. 
The results show that mole fraction 0 2 and He-based jet mixture fractions correlate positively 
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but exhibit significant deviations from the mass fraction C-based localized jet mixture fraction. 
In general, the 0 2 -based mixture fraction underpredicts the C-based mixture fraction due to jet 
dilution and combustion, with additional discrepancies near § = 1 due to non-equilibrium 
conditions in the RQL quick mixing zone. The He tracer, on the other hand, overpredicts the 
C-based mixture fraction, possibly due to differences in density and diffusivity, and exhibits 
significant scatter most likely attributable to differences in gas density and turbulent diffusivity. 
However, the data show a much better quantitative agreement between the 0 2 , He, and C-based 
methods when assessing planar mixing fields in terms of spatial unmixedness. Although the 
combustor rig was designed to test RQL combustion, the results of this experiment can 
potentially be applied more generally as a diagnostic to assess air- fuel mixing in other types 
of reacting systems. 
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TABLE 1.— FORMULAE AND ASSUMPTIONS FOR CALCULATING JET MIXTURE FRACTION USING 

DIFFERENT CONSERVED SCALARS 


Conserve 
d Scalar 

Jet mixture fraction (f) 
formula 

Known 

quantities 

Unknown 

quantities 

Equations and assumptions 

ic r G 

ysample 

f =1 lc 

J Y C y crossflow 

y crossflow ysample 

r _ 1 0 ~ 1 0 

J Yq y crossflow y jet 

x o ~ 1 o 

X CO, dry, X C 02,dry, 
X 02, dry, X TH C,dry 

X CO, wet, Xc02,wet , 

Xo2,wet , Xc3H8,wet > 
Xc2H4,wet , X H 20,wet, 
Xm,wet,Xm,wet 

1-4) X i dry = , i = CO, C0 2 , 0 2 , THC 

1 — X' 11,0 

5) C:H = 3:8 

6) 0:N = 0.209/0.791 

7) X H2 = 0.65A CO 

8) Z^=1 

i 

9) THC = C 3 H 8 + C 2 H 4 

Y02 

ysample 
f _ °2 

y o 2 yjet 

1 o 2 

X CO, dry, Xc02,dry , 
X 02, dry, X THCt dry 

X CO, wet, Xc02,wet , 
X 02, wet, Xc3H8,wet > 
Xc2H4,wet, X H 20,wet, 
Xjs[2, wet,Xn2, wet 

1 -9) same as above 

10) xy fl0 " = 0 ; Xg = 20.9% 

Y„ e 

ysample 

f =1 He 

J Y He y crossflow 

X CO, dry, X C 02,dry, 
X 02, dry, , X TH Qdry, 
XHe, dry 

X CO, wet, Xc02,wet , 

X 02, wet, Xc3H8,wet , 
Xc2H4,wet , X H 20,wet, 
X]Sf2,wet , Xh2, wet, Xn e , wet 

1-9) same as above 

10) xg=o 

X 02 

ysample \(J 

p Q2 L z A sam pi e 

Jx °, - X je< - 20.9% 

X 02, dry 

None 

1) X c ™ sfl ™ = 0 ; Xg = 20.9% 

2) M sample ~ constant 

2) X (>2 iry oc X ()2 wel 

x„ e 

ysample 

f =1 He 

J X He crossflow 

[He] 

| L J sample 

\^ e \Max{ x / R =- 1 ) 

XHe, dry 

None 

1) X£=0 

2) Plane x/R = -1 is approximately uniform) 

3) M samp i e ~ constant, 

4 ) X He, dry ^ ^ He, wet 



TABLE 2.— KEY PROPERTIES OF SPECIES USED IN THE 
DETERMINATION OF JET MIXTURE FRACTION 


Compound 

Molar mass, 
(g/mol) 

Maximum 

concentration, 

(%, dry) 

Maximum 

density, 3 

(g/m 3 ) 

Air 

28.8 

100.0 

931.7 

CO 

12.0 

10.0 

151.2 

co 2 

44.0 

13.2 

180.0 

0 2 

32.0 

20.9 

273.5 

C 3 H 8 

44.1 

2.2 b 

38.8 

h 2 o 

34.0 

14.8 

205.2 

He 

4.0 

0.3 

0.5 

Ne 

20.2 

0.2 

1.9 

Ar 

39.9 

— 

— 


PM 

Density (/?) is given by p = 100 C , where P = 1 atm, T = 298 K, 

RT 

M = molar mass, R = 8.3 14 J/molK, and C = concentration in %. 
b Measured as total unburned hydrocarbons corrected to propane 
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TABLE 3.— OVERALL CONSTITUENT MASS FRACTIONS IN 
ZONES UPSTREAM AND DOWNSTREAM OF JET MIXING 
SECTION, BASED ON FLOW RATES OF FUEL, AIR, 

AND TRACER GAS 


Zone 

Helium 

Propane 

Air 

y a 

1 rich 

4.4 x 10' 4 

0.097 

0.90 

y b 
1 lean 

1.3 x 10' 4 

0.028 

0.97 

a rich = helium + propane + air 
b lean = rich + air jets 


TABLE 4.— STANDARD ERROR (a) AND COEFFICIENTS OF 
DETERMINATION (R 2 ) FOR 0 2 AND He-BASED MIXTURE 
FRACTIONS (f x ) AND SPATIAL UNMIXEDNESS (U s ) 
WITH RESPECT TO f Y c AND U s , Y c 


Parameter 

fxOl 

fxHe 

Us,X02 

Us,XHe 

a 

0.09 

0.12 

0.06 

0.08 

R 2 

0.94 

0.90 

0.96 

0.93 


a The a and R 2 values for jet mixture fraction / were calculated using all the 
data obtained for each module, except for points at the upstream plane x/R = -1 
where /is uniform and approximately zero. 
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Overall lean products 


x/R = 



2R/3 N 
R/3 S 12\ 

N 7 13 'n 

2 8 

1 4 g 14 1 ^Orifices 

5 


6 


10 

\ 15 

16/ 


Rich reacting flow 

Figure 2. — Location of planes and sampling points in RQL combustor module. 
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Figure 3. — Setup of helium injection system, sampling train, and analysis system. 



Figure 4. — Comparison between carbon mass fraction-based (fyc) and 
oxygen concentration-based (fxo2) i et mixture fractions. 
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Figure 5. — Comparison between carbon mass fraction-based (fyc) and 
helium concentration-based (fxHe) i et mixture fractions. 
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Figure 6. — Carbon mass fraction-based (fyc) mixture fraction fields. 
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Figure 7. — Oxygen concentration-based ( fx02 ) mixture fraction fields. 
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Figure 8. — Helium concentration-based (fxHe) mixture fraction fields. 
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Figure 9. — Comparison between carbon mass fraction-based (Us,yc)> 
oxygen concentration-based (Us,x02)> ancl helium-concentration 
based ( Us,X,He ) spatial unmixedness for all modules. 
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Figure 10. — Spatial unmixedness (Us) values at each measurement plane. 
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